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CONVERGENCE  PERFORMANCE  OF  ADAPTIVE  DETECTORS 

PART  3 


1.0  INTRODUCTION 

In  Ref.  1,  two  schemes  for  adaptive  detection  were  compared:  Kelly’s  generalized  likelihood 
ratio  test  (GLRT)  [2]  and  the  mean  level  adaptive  detector  (MLAD).  Detection  performance  Pp  was 
predicted  for  the  two  schemes  under  the  assumptions  that  the  input  noises  are  zero-mean  Gaussian 
random  variables  that  are  temporally  independent  but  spatially  correlated,  and  the  desired  signal's 
amplitude  is  Rayleigh  distributed.  Pp  was  computed  as  a  function  of  the  false  alarm  probability,  the 
number  of  input  channels,  the  number  of  independent  samples-per-channel,  and  the  matched  filtered 
output  signal-to-noise  (S/N)  power  ratio.  The  GLRT  was  shown  to  have  better  detection  performance 
than  the  MLAD.  The  difference  in  detection  performance  increased  as  one  used  fewer  input  samples; 
however,  the  required  number  of  samples  necessary  to  have  only  a  3  dB  detection  loss  for  both  detec¬ 
tion  schemes  is  approximately  the  same.  This  is  significant  since,  for  the  present,  the  MLAD  is  con¬ 
siderably  less  complex  to  implement  than  the  GLRT. 

The  general  problem  of  signal  detection  in  a  background  of  Gaussian  noise  for  an  adaptive  array 
was  first  addressed  by  Kelly  [2]  by  using  the  techniques  of  statistical  hypothesis  testing.  In  Ref.  2, 
the  problem  is  formulated  as  a  binary  hypothesis  test  where  one  hypothesis  is  noise  only  and  the  other 
is  signal-plus-noise.  A  given  input  data  vector  (called  the  primary  data  vector)  is  tested  for  signal 
presence.  Another  set  of  signal-free  data  vectors  (called  the  secondary  data  vectors)  is  available  that 
shares  the  unknown  covariance  matrix  M  of  the  noise  in  the  primary  data  vector.  A  likelihood  ratio 
decision  rule  was  derived,  and  its  performance  was  evaluated  for  the  two  hypothesis. 

Kelly’s  detector  uses  the  maximum  likelihood  (ML)  estimates  for  the  unknown  parameters  of  the 
likelihood  ratio  test  (LRT).  The  unknown  parameters  are  the  spatial  covariance  matrix  and  the 
unknown  signal’s  complex  amplitude  (assumed  in  Kelly’s  analysis  to  be  a  nonrandom  constant).  This 
detection  scheme  is  commonly  referred  to  as  the  GLRT  and  is  referenced  in  this  report  as  such. 

A  less  complex  adaptive  detection  scheme  is  found  by  implementing  MLAD.  The  MLAD  is 
essentially  an  adaptive  matched  filter  (AMF)  followed  by  a  mean  level  detector  (MLD)  13.4).  Input 
samples  used  in  determining  the  MLD  threshold  are  derived  from  a  block  of  data  passing  through  the 
AMF.  This  same  block  of  data  is  used  to  calculate  the  AMF  weights.  The  squared  magnitude  of 
each  of  these  same  samples  as  processed  through  the  AMF  is  u^ed  as  a  test  statistic  and  compared 
against  an  MLD  threshold  (an  average  of  the  instantaneous  powers)  that  does  not  contain  the  given 
test  statistic  sample.  We  further  clarify  the  implementation  terminology  by  calling  this  an  MLAD 
with  concurrent  data  samples.  In  Ref.  5,  an  analysis  was  performed  for  an  MLAD  with  noncon¬ 
current  data,  i.e.,  the  block  of  data  that  passes  through  the  AMF  that  is  used  to  determine  the  MLD 
threshold  is  statistically  independent  of  the  block  of  data  used  to  calculate  the  AMF  weights. 

We  note  for  both  Kelly’s  GLRT  and  the  MLAD  that,  under  the  above  stated  assumptions,  the 
PF  does  not  depend  on  M  (a  second  order  characterization  of  the  external  noise  environment).  Hence 
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detectors  exhibit  the  desirable  constant  false  alarm  rate  (CFAR)  property  of  having  the  P y  be  indepen¬ 
dent  of  the  covariance  matrix. 

Here  we  extend  the  results  of  Ref.  1  to  include  the  effects  of  what  we  term  "desired  signal  con¬ 
tamination.”  By  this,  we  mean  that  a  significant  level  of  a  desired  signal  is  present  in  the  secondary 
data  vectors,  which,  for  both  Kelly's  GLRT  and  the  MLAD.  are  used  to  estimate  the  unknown 
covariance  matrix  M.  We  assume  that  the  contaminating  desired  signal  is  statistically  independent  of 
the  desired  signal  in  the  primary  data  vector.  In  practice,  contamination  can  be  caused  by  a  variety  of 
mechanisms,  such  as  for  the  radar  problem,  multiple  target  returns  being  present  at  different  ranges  of 
the  pulse  compressor  output.  In  this  case,  a  number  of  the  multiple  target  returns  at  distinct  times 
that  have  the  desired  signal’s  waveform  may  be  in  the  secondary  data. 

For  our  analysis,  we  simplify  the  contamination  model  by  only  considering  a  single  source  of 
contamination.  In  addition,  this  single  source  contaminates  only  one  sample  vector  of  the  secondary 
data  vectors.  Two  performance  measures  are  affected:  Pp  and  Py.  Because  Ph  is  normally  chosen 
under  the  condition  of  no  contamination,  we  redesignate  the  probability  of  false  alarm  in  the  presense 
of  contamination  as  the  ghosting  probability  PG.  We  do  this  since  ghosi>  (in  the  radar  sense)  are 
desired  signal-induced  false  detections. 

Here  we  present  results  on  the  detection  and  ghosting  performance  of  GLRT  and  MLAD  in  the 
presence  of  contamination.  As  in  Ref.  1.  we  assume  that  the  complex  desired  signal  amplitude  is  a 
complex  zero-mean  Gaussian  random  variable  (r.v.)  of  unknown  variance  with  independent  and  iden- 
tcally  distributed  (i.i.d.)  real  and  imaginary  parts  (the  magnitude  of  this  amplitude  is  Rayleigh  dis¬ 
tributed).  Under  the  GLRT,  we  would  have  to  reformulate  Kelly's  detector  with  the  variance  of  the 
unknown  signal  amplitude  as  an  unknown  parameter,  and  find  the  ML  estimate  of  this  quantity.  This 
proved  to  be  mathematically  tedious.  In  lieu  of  implementing  this  new  GLRT.  we  choose  to  evaluate 
Kelly’s  GLRT,  as  it  is  defined  in  his  paper.  As  noted  by  Kelly,  no  optimality  properties  are  claimed 
for  this  test.  The  form  of  the  test  is.  however,  reasonable. 

2.0.  GENERALIZED  LIKELIHOOD  RATIO  TEST 

2.1  Detector  Form 

A  mathematical  formulation  of  the  adaptive  detection  problem  that  leads  to  the  GLRT  is  given 
by  Kelly  (2).  We  now  summarize  that  formulation.  Two  sets  of  input  data  are  used,  called  the  pri¬ 
mary  and  secondary  inputs.  The  secondary  inputs  are  assumed  not  to  contain  the  desired  signal.  Set 

X  =  N  X  K  matrix  of  secondary  input  data.  The  nth  row  represents  the  K  samples  of  data 

on  the  nth  channel,  where  n  =  1.2 . N.  The  samples  in  the  Ath  column  are 

assumed  time-coincident. 

x  =  primary  data  vector  of  length  N. 

s  =  desired  steering  vector  of  length  N. 

Consider  the  two  hypothesis: 


A/,,  :  x  =  n.  and  <  1  > 

:  x  =  n  +  (/  s.  (2) 
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where  H0  is  the  noise  only  hypothesis,  n  is  a  noise  vector  of  length  N,  H  {  is  the  signal-plus-noise 
hypothesis,  and  a  is  the  unknown  complex  signal  amplitude.  We  make  the  following  assumptions: 

(Al)  Input  noises  are  complex  zero-mean  stationary  Gaussian  r.v.'s.  The  real  and  imaginary 
parts  of  a  given  input  noise  sample  are  i.i.d.  with  respect  to  each  other  (an  r.v.  with  these 
characteristics  is  called  a  circular  Gaussian  process). 

(A2)  Input  noise  samples  are  temporally  statistically  independent. 

(A3)  The  secondary  data  is  statistically  independent  of  the  primary  data. 

(A4)  K  >  N. 

The  GLRT  is  formulated  as  follows.  The  joint  probability  density  function  (PDF)  under  each 
hypothesis  over  all  measured  data  is  found.  For  this  problem,  this  is  straightforward,  since  the  sam¬ 
ple  vectors  are  assumed  independent  and  each  vector  has  an  associated  (V-dimensional  Gaussian  PDF. 
If  there  are  any  unknown  parameters,  the  PDF  of  the  inputs  is  maximized  over  all  unknown  parame¬ 
ters  separately  for  each  of  the  two  hypotheses.  The  maximizing  parameter  values  are  by  definition 
the  ML  estimators  of  the  parameters.  Hence  the  maximized  PDFs  are  obtained  by  replacing  the 
unknown  parameters  by  their  ML  estimates.  The  ratio  of  the  resultant  maximum  of  PDFs  is  found 
(the  ratio  of  the  PDF  under  //[  to  the  PDF  under  H0).  This  ratio  is  checked  to  see  if  it  exceeds  a 
preassigned  threshold  t. 

Kelly  shows  that  the  GLRT  for  the  adaptive  detection  problem  is  given  by 

H\ 

| sH  Rx  *  x | 2  > 

- - ; - ^  t,  (3) 

(sH  Rx  s)(l  +  xH  Rx  x]  < 

H0 


where 


Rx  =  XXH,  (4) 

and  H  denotes  the  conjugate  transpose  matrix  operation.  We  recognize  Rx  as  proportional  to  the  ML 
estimate  of  the  input  covariance  matrix.  We  note  also  that  the  desired  signal’s  unknown  complex 
amplitude  a  has  been  estimated  and  is  accounted  for  in  Eq.  (3).  The  elements  of  Rx  are  r.v.'s  that 
are  functions  of  the  input  samples  (the  elements  of  X).  It  is  straight-forward  to  show  that  the  proba¬ 
bilistic  measure  of  the  set  of  X ,  for  which  Rx  is  singular,  is  zero.  Hence,  when  assessing  detection 
performance,  we  can  always  assume  that  Rx  (or  any  other  matrix  that  has  the  form  given  by  Eq.  (4)) 
is  invertible. 

For  the  signal  contamination  model,  we  make  the  additional  assumption: 

(A5)  A  statistically  independent  desired  signal  is  always  present  in  the  secondary  data  under  H{) 
or  H  | .  It  is  only  present  on  the  Kth  sample  vector  (the  Kth  column  of  X). 

For  (A5),  the  choice  of  which  secondary  sample  vector  in  time  is  contaminated  is  arbitrary  and  does 
not  affect  the  analysis. 
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2.2  Statistically  Equivalent  GLRT 

In  Ref.  1,  a  statistically  equivalent  GLRT  was  derived  that  was  used  to  formulate  in  simple 
fashion  the  Pp  and  PF  probabilities  of  the  adaptive  detector  with  no  signal  contamination.  With  one 
modification  of  the  development  in  Ref.  1,  we  can  derive  a  statistically  equivalent  GLRT  that  will 
allow  us  to  formulate  the  Pp  and  PG  for  the  adaptive  detector  with  signal  contamination  (which  is 
modeled  as  given  by  (A5)).  We  briefly  outline  the  methodology  used  to  obtain  the  statistically 
equivalent  GLRT. 

As  in  Kelly’s  (and  Reed,  Mallet,  and  Brennan  16])  development,  we  can  matrix  transform  the 
input  vectors  by  an  N  x  N  matrix  A,  which  has  the  properties  that  the  input  noise  vectors  are  spa¬ 
tially  whitened,  each  input  element  has  noise  power  normalized  to  one,  and 

As  =  (0,0,  •0,(sH  M s)l/2)  s  so,  (5) 

where  all  of  the  desired  signal  has  been  placed  into  N  th  channel  (note  that  in  Kelly’s  paper,  the  signal 
was  placed  into  the  first  channel;  for  our  analysis,  we  place  the  signal  into  the  N  th  channel). 

In  addition,  set 

z  =  Ax,  and  (6(a)) 
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where 


(10) 


and  «u,«2i 'm22’vi ,  and  v’2  are  independent  r.v.’s.  Furthermore,  un  has  the  \  PDF  of  order 
2 (K  -  N)  with  <r  =  0.5,  v,  has  the  x  PDF  of  order  2 (N  -  1)  with  <r  =  0.5  under  H , ,  v2  is  the 
sum  of  the  desired  signal  in  the  primary  data  plus  a  circular  Gaussian  r.v.  with  power  equal  1,  and 
u 2i  is  a  circular  Gaussian  r.v.  with  power  equal  1.  In  addition,  if  we  adopt  the  notation  of  Ref.  1. 
then 


U22 


k  =,V 


(11) 


where  cJw  k  =  N,N  +  1  -  1  are  i.i.d.  circular  Gaussian  r.v.  with  power  equal  1  and 

cJva.'-  1  ’  is  the  sum  of  the  contaminating  desired  signal  plus  a  circular  Gaussian  r.v.  with  power  equal 

1. 

We  use  the  statistically  equivalent  GLRT  given  by  Eq.  (9)  and  the  aforementioned  probabilistic 
characterizations  of  u  M  ,u2\  .«22»v  i  a°d  v2  as  the  starting  point  for  our  analysis. 

2.3  Probability  of  Detection 

Under  the  H\  hypothesis,  we  assume  that  the  primary  vector’s  desired  signal’s  amplitude  (or 
magnitude)  is  Rayleigh  distributed  and  the  signal’s  phase  is  uniformly  distributed  between  (0,2ir). 
This  implies  that  the  desired  signal  itself  is  a  complex  circular  Gaussian  r.v.  ^Let  the  desired  signal’s 
input  power-per-channel  before  any  matrix  transformation  be  equal  to  o~ .  After  the  A  matrix 
transformation  (whitening,  normalizing,  and  placing  the  signal  into  the  N  th  channel),  the  signal 
power  in  the  Ath  channel  is  a;  =  o~  sH  M-Is.  Thus  under  H\,v2  is  a  complex  circular  Gaussian 
r.v.  with  power  equal  to  +  1,  where  the  1  represents  the  noise  power-lcvel-per-ehannel  after  the  A 
matrix  transformation.  In  similar  fashion,  let  the  contaminating  desired  signal’s  input  power-per- 
channel  before  any  matrix  transformation  be  equal  to  w  Again  we  assume  the  contaminating  desired 
signal’s  amplitude  is  Rayleigh  distributed  with  phase  uniformly  distributed  between  (0. 27t).  After  the 
A  transformation,  a *  =  a,  sH  M-1  s.  Thus  under  H{)  or  H j ,  z{\k  1 1  is  a  complex  circular  Gaussian 
r.v.  with  power  equal  to  a~  +  1,  where  again  the  1  represents  the  noise  power-level-per-channel  after 
the  A  matrix  transformation. 


Recall  that  u2\  is  a  complex  circular  Gaussian  r.v.  with  power  equal  to  1.  We  can  rewrite  Eq. 

(9)  as 

Hx 

M2  <  T\  ,12) 

Ho 

where 
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It  is  straightforward  to  show  that  a,  when  conditioned  on  uM  and  Vj,  is  a  complex  circular  Gaussian 
r.v.  with  power  equal  to  1.  It  is  well  known  [7]  that  the  conditional  probability  of  detection  is  given 
by 


P(D  |  uu,u22,Vi)  =  e  T'  =  exp 


«11  +  V2\ 
u2u(o2  +  1)  + 


05) 


where  P{D  \  ,  )  denotes  the  conditional  probability  of  detection. 
We  set 


V  =  u  22’ 

(16) 

M  =  «11. 

(17) 

v  =  vj,  and 

(18) 

V 

(19) 

Then  Eq.  (15)  becomes 


P(D  |  r,ri)  =  exp  <  - 


1  +  r 


<jj  +  1  +  r 


VT 


(20) 


The  PDFs  of  n  and  N  are  x2  of  order  2 (K  —  N  +  2)  and  2 (N  -  1),  respectively,  with  o2  =  0.5  and 
are  given  by 


pM  = 


l 


(K  -  N  +  1)! 


HK  n  +  i  e  M  yfl  >o,  and 


p..(v)  =  - - -  vN  ~  2  e~v  ,v  >  0. 

(N  -  2)! 

The  PDF  of  rj  is  derived  in  the  Appendix  and  is  given  by 

/-,<■»  =  — i - ' - vK~"e  l' 

’  (a2.  +  1  )(K  ~  N  -  1)!  Jo 


(21) 


(22) 


L  x*  N~'  exp  i  - — 


a r. 


ar.  +  1 


>d\ 


.r]  >  0. 

By  using  elementary  probability  theory,  it  is  straightforward  to  show  that 

0 


(23) 


Prir)  =  L  -4-  PvW)p„  \— 

0  l r  J 


<10. 


(24) 
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By  inserting  expressions  for  pM  and  pv  as  given  by  Eq.  (21)  and  (22),  respectively,  and  simplifying 
results  in  the  following  expression  for  pr: 


Pr(r) 


K\ 


,N-  2 


(N  -  2 )\(K-N  +  1)!  (l  +  r) 


K  +  1 


,r  >  0. 


(25) 


If  we  set  q  =  1/(1  +  r),  it  is  straightforward  to  show  that 


Pq(<i)  = 


K\ 


(M  —y\ukr  —  N 


I  1  \  I 


(1  -q)N~2qK-N  +  ' 


(26) 


which  is  the  PDF  derived  by  Reed  et  al.  (6]  for  the  normalized  instantaneous  S/N  power  ratio  and 
which  results  if  the  sampled  matrix  inversion  (SMI)  algorithm  is  used.  By  substituting  q  for  r  in 
Eq.  (20), 


p(D\q,rj)  =  exp  - 


1 

q  a2  +  \ 


VT. 


(27) 


If  we  multiply  Eq.  (27)  by  the  PDF  of  tj  as  given  by  Eq.  (23)  and  integrate  over  the  support  of  r /,  it 
is  straightforward  to  show 


Thus 


p{D  | q)  = 


K-N 

ol  +  1 


\K~N-'d\ 

T  1  +  \a2c 

qa2s  +  l  a2.  +  1 


K-N  +  1 


(28) 


P(D)  =  P(D  |  q)  pq(q)dq ,  (29) 

where  pq(q)  is  given  by  Eq.  (26).  We  set  (S/N)  op,  —  a2,  where  (S/N)„p,  is  the  optimal  S/N  output 
power  ratio  of  the  matched  Filter  ( K  =  <x>).  We  can  write  this  in  this  way  because  the  output  noise 
power  of  the  TV th  channel  has  been  normalized  to  l  and  the  output  of  the  Nth  channel  is  the  optimal 
matched  filter  output.  We  also  set  (S/N)a,„  =  a2,  where  (S/N)(W1  is  the  contaminated  S/N  output 
power  ratio  of  the  matched  filter  ( K  =  oo). 


2.4  Probability  of  Ghosting 

The  probability  of  a  ghost  is  easily  derived  from  Eqs.  (28)  and  (29)  by  setting  a;  -  0.  It  is 
found  that 


Pc 


K-N 

<7^  +  1 


\K~N-'d\ 

1  +  \a2/  *~'V  +  I 
1  +  a2 


If  we  set  a2  =  0,  then 


Pg  =  Pf 


1 

( T  +  l)*~"  +  l  ’ 


(30) 


(31) 
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which  is  the  probability  of  false  alarm  for  the  GLRT  derived  in  Ref.  I  under  the  condition  that  there 
is  no  contamination. 

3.0  MEAN  LEVEL  ADAPTIVE  DETECTOR 

3.1  Detector  Form  and  Statistically  Equivalent  MLAD 

In  Ref.  1,  an  MLAD  was  formulated  for  the  adaptive  detection  problem.  In  this  subsection,  we 
briefly  review  this  formulation  for  the  MLAD  and  give  a  statistically  equivalent  form  that  will  be 
used  to  obtain  results  for  PD  and  PG. 

An  intuitive  form  of  adaptive  detection  is  found  by  implementing  the  MLAD.  The  MLAD  is 
essentially  an  AMF  followed  by  an  MLAD,  as  illustrated  in  Fig.  1.  The  MLAD  is  designed  to  per¬ 
form  detections  over  a  block  of  data  by  using  just  this  block  of  data  in  determining  the  AMF  weights 
and  the  MLD  threshold.  The  MLAD  works  as  follows.  Let  there  be  N  channel  and  K  +  1  samples 
per  channel.  Define 

x  =  primary  N  -length  data  vector; 

=  secondary  N  -length  data  vector,  k  =  1.2, . . .  ,K , 

XaUg  =  (x  |  X)  =  augmented  N  x  (K  +  1)  matrix  of  input  data,  and 

^0  =  -^aug  ^aug- 

The  N  -length  weighting  vector  w  for  the  AMF  is  found  by  using  the  SMI  algorithm  and  is  given  by 

w  =  R0  1  s.  (32) 

This  weight  is  used  in  the  detection  rule  given  by 


To  £  |wHx*)2,  (33) 

k  =  I 


probability.  We  see  that  Eq.  (32)  is  the  algorithmic 
representation  of  the  AMF  and  Eq.  (33),  the  MLD. 

Note  that  we  have  included  the  primary  data  vector  in  the  R0  estimate  and,  hence,  in  the  w  esti¬ 
mate.  In  a  practical  situation,  this  might  be  done,  since  it  is  more  numerically  efficient  to  compute 
one  weighting  vector  over  the  entire  data  block  than  it  is  to  compute  a  distinct  weighting  vector  for 
each  point  in  the  block.  However,  the  presence  of  the  desired  signal  (under  H !)  will  affect  detection. 
In  Eq.  (33),  the  primary  data  vector  is  varied  across  the  K  +  1  data  snapshots,  where  the  x^  used  on 
the  right  side  of  Eq.  (33)  does  not  include  a  selected  primary  data  vector. 


H 


"H  |  2  > 
w  x|z  < 


where  To  is  chosen  to  control  the  false  alarm 
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Furthermore,  it  is  shown  in  Ref.  1  that  a  statistically  equivalent  MLAD  is  given  by 


M  1 1  v2  ~  'T«21 


2  2  t 

U  |  |  U  n  I  \  , 


where  u  n  ,u2\ .m22  AM ,  and  v2  are  as  probabilistically  characterized  in  Section  2.2. 
We  note  for  T,  >0  that  if 


then  Hn  is  declared. 


3.2  Probability  of  Detection 


Again  under  the  H ]  hypothesis,  we  assume  the  primary  vector's  desired  signal  and  the  contam¬ 
inating  desired  signal  are  complex  circular  Gaussian  r.v.  (the  amplitude  is  Rayleigh  distributed).  As 
in  our  analysis  of  the  GLRT,  a)  =  b~  sH  M“‘s  and  a2.  =  a~  sw  M  1  s.  Assume  Eq.  (38)  is  not 
true.  We  can  write  the  decision  rule  given  by  Eq.  (37)  as 


2  >r 
<  7> 


where 


M  1 1  v  2  ~  v  1  **  21 

(«ll  (<*?  +  1)  +  v])U1  ' 


T[  =  - 5 ; - r]T , , 

{q  -  Tx)(qo]  +  1) 


r)  is  defined  by  Eq.  (16),  and  q  =  (1  +  r)  *.  As  before,  under  H  j ,  v  j  is  a  complex  circular 
Gaussian  r.v.  with  power  equal  to  a]  +  1,  and  «2|  is  the  same  with  power  equal  to  1.  Furthermore, 
a  is  the  same  with  power  equal  to  1 .  Thus 


P(D\q,V) 


0  if 


q  <  T\ 
otherwise. 


The  PDFs  of  r)  and  q  are  given  by  Eqs.  (23)  and  (26).  respectively.  If  we  multiply  P\D  !</.>/)  by 
p,}(r\)  and  integrate  over  the  support  of  r 7.  after  some  simplification,  it  will  be  found  that  for  q  >  7  , . 
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P(D  |  4) 


K-N 
or.  +  1 


!„ 


.  K  -  \  \ 


d\ 


1  +  \o~ 


+ 


(q  -  T ,  )(qo~  +  1)  \  +  a~ 


K  -  .V  +  I 


(43) 


Because  P(D  |  q)  =  0  for  q  <  T{ ,  it  follows  that 

PD  =  j t,  P(-D  I  <?)  (i<1'  <44> 

-> 

where  pq(q )  is  given  by  Eq.  (26).  Again  we  set  tS/N)<y),  =  a;  and  (S/N),.,,,,  =  a~ . 

3.3  Probability  of  Ghosting 

The  PF  is  found  by  setting  (S/N)()/„  =  0  in  Eq.  (43).  The  following  equations  result: 


,k  -s  - 1 


P(G\q)  - 


K-N  f  i 
ct?  +  1  0 


d\ 


+ 


I  4-  X(7" 
1  +  or. 


K -S  +  l 


.  q  >  T , 


(45) 


.q  ^  t i 


and 


PG  =  \T  P(G  |  q)pc/(q)dq. 


If  we  set  a =  0,  it  is  straightforward  to  show  that 

Pc  =  Pf  =  (1  "  T\)K. 


(46) 


(47) 


which  is  the  false  alarm  probability  of  the  MLAD  derived  in  Ref.  1  under  the  condition  that  there  is 
no  contamination. 


4.0  RESULTS 

Here  we  present  results  for  tr..1  detection  probability  in  contamination  PD  and  ghosting  probabil¬ 
ity  PG  for  the  GLRT  and  MLAD  vs  the  independent  parameters:  the  quiescent  prob¬ 
ability  of  false  alarm  (when  there  is  no  contamination)  PF:  the  steady  state  (K  =  oo)  s/N  output 
power  ratio  of  the  matched  filter  (S/N),Y„:  the  contaminated  S/N  power  ratio  (S/N)(.„„;  the  number  of 
independent  samples  per  channel  K  of  secondary  data;  and  the  number  of  input  channels  N.  We  set 
K  =  MN,  where  M  is  a  positive  integer  >  2.  and  use  M  instead  of  K  as  an  independent  parameter 
called  the  degrees-of-freedom  factor. 

Because  there  are  many  independent  parameters,  we  present  results  for  a  representative  set  as 
shown  in  Figs.  2  through  19.  Here  we  plot  PD  and  PG  vs  (S/N),„„  for  N  =  2,  5.  10.  30; 

PF  =  10  6,  10~10,  and  M  =  2,  3 . 6.  For  plotting  Pn.  we  choose  (S/N)(l/„  to  equal  either  20  dB 

or  30  dB,  where  a  given  (S/N),y„  will  yield  a  steady  state  (K  =  oo)  Pn,  which  is  indicated  by  the 
horizontal  line  in  Figs.  10  through  19.  For  plotting  PG  in  the  steady  state.  PG  =  Pf .  which  is  indi¬ 
cated  by  horizontal  line  in  Figs.  2  through  9. 


I! 


S/N  CONTAMINATED  (dB) 

3  —  Kelly  detector:  PG  tor  contaminated  signal  ,V  -  5.  Ph 


I  7)  6 


PROBABILITY  OF  GHOSTING  (PG) 
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Fig.  4  —  Kelly  detector:  PG  for  contaminated  signal  N  =  10.  PF  =  l.D-6 


Fig.  5  —  Kelly  detector:  PG  for  contaminated  signal  /V  =  30.  PF  -  I  D-b 
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Fig.  16  —  Mean  level  detector:  PD  contaminated  N  =  2.  PF  =  1.0-6.  S/N  opt  =  20  dB 


opt  =  30  dB 


S/N  CONTAMINATED  (dB) 

Mean  level  detector:  PD  contaminated  N  =  5.  PF  =  1.0-6.  S/N 


PROBABILITY  OF  DETECTION  WITH  CON.  (PDC) 
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We  found  that  the  PG  performance  results  were  identical  for  both  the  GLRT  and  MLAD  and  a 
given  set  of  input  parameters.  Thus,  we  give  only  the  PG  performance  results  for  the  GLRT.  Some 
pertinent  observations  to  be  made  from  Figs.  2  through  19  are: 

1.  Both  PD  and  PG  degrade  monotonically  with  increasing  contamination  (S/N),.,,,,. 

2.  A  small  amount  of  contamination  ((S/N),.„„  =  -  10  dB)  decreases  PG.  However.  10  dB  of 
contamination  is  necessary  make  PG  decrease  approximately  by  a  factor  of  10. 

3.  For  most  cases,  Pd  begins  to  significantly  decrease  in  10  to  20  dB  of  contamination. 

4.  For  small  M  (2-6),  the  PG' s  are  within  a  factor  of  five.  This  spread  increases  with 
increasing  N. 

5.  As  noted  in  Ref.  1,  PD  monotonically  increases  with  increasing  M.  For  (S/N),.,,,,  =  0.  the 
steady  state  Pp' s  are  indicated  by  the  flat  region  of  each  curve. 

6.  For  Pd.  the  GLRT  and  MLAD  have  similar  relative  performance  trends. 

One  of  the  more  significant  results  indicated  by  the  curves  is  that  the  ghosting  probability  does 
not  increase  in  the  presence  of  contamination.  Hence,  the  CFAR  capability  of  both  the  GLRT  and 
MLAD  is  not  degraded  in  the  sense  that  the  false  alarm  probability  (with  contamination)  is  upper 
bounded  by  the  quiescent  false  alarm  probability  (no  contamination).  However,  because  PG  is 
decreasing,  inherently  the  detector’s  (MLAD  or  GLRT)  variable  threshold  is  increasing  and,  hence, 
the  detection  probability  decreases. 

One  final  note.  A  simple  solution  to  decrease  the  effects  of  signal  contamination  is  to  use  a 
large  number  of  samples  (make  K  or  M  large).  However,  we  caution  against  this  solution  in  that  this 
was  a  simplified  analysis  where  only  one  sample  vector  was  contaminated.  Obviously,  if  we  take 
enough  samples,  the  effects  of  this  one  sample  can  be  significantly  diminished  by  averaging  over 
many  samples.  In  reality,  making  K  large  can  result  in  even  more  signal  contamination,  since  there 
may  be  more  opportunities  for  this  to  occur. 

5.0  SUMMARY 

Two  schemes  for  adaptive  detection,  Kelly’s  generalized  likelihood  ratio  test  (GLRT)  and  the 
mean  level  adaptive  detector  (MLAD),  have  been  analyzed  with  respect  to  the  deleterious  effect  of 
desired  signal  contamination  of  the  data  used  to  compute  the  sampled  covariance  matrix  for  the  two 
detections.  This  effect  can  occur  when  more  than  one  desired  signal  is  present  in  the  sampled  data. 
Detection  probability  Pd  and  false  alarm  performance  (ghosting  probability  PG)  were  predicted  for 
the  two  schemes  under  the  assumptions  that  the  input  noises  were  Gaussian  random  variables  that 
were  temporally  independent  but  spatially  correlated;  and  the  desired  signal's  amplitude  was  Rayleigh 
distributed.  Pp  and  PG  were  computed  as  a  function  of  the  false  alarm  probability  with  no  contami¬ 
nation  PF.  the  number  of  input  channels,  the  number  of  independent  samples-per-channel.  the 
matched  filtered  output  S/N  power  ratio,  and  the  S/N  of  the  contaminating  desired  signal.  The  Pp 
and  Py  were  obtained  for  a  number  of  representative  cases. 

It  was  found  that  both  Pp  and  PG  decreased  with  increasing  levels  of  contamination.  The  P(, 
performance  was  almost  identical  for  the  GLRT  and  MLAD.  The  Pp  performance  for  the  two  adap¬ 
tive  detectors  showed  similar  relative  performance  trends.  Significantly,  it  was  shown  that  the  ghost¬ 
ing  probability  does  not  exceed  Py  in  the  presence  of  contamination.  Hence,  the  CFAR  capability  of 
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the  GLRT  or  MLAD  is  not  degraded  in  the  sense  that  false  alarm  probability  is  upper-bounded  by  the 
quiescent  false  alarm  probability  (no  contamination). 
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Appendix 

THE  PROBABILITY  DENSITY  FUNCTION  (PDF)  OF  r, 


Starting  with  Eq.  (11),  we  write 


=  =  £  UJfi-'M2  +  I4VU|2 


where  k  =  N  + 1,  .  .  .,k  - 1  are  independently  and  identically  distributed  circular  Gaussian 

random  variable  (r.v.)  with  power  equal  to  1,  and  cl,*  "  is  a  circular  Gaussian  r.v.  with  power 
equal  to  a"  +  1 .  Define 


-v  -  I  Z$K 


i  - 


The  PDFs  of  .t  and  y  are  given  by 


K-  1 

|  ,(A' -  I)  |  2 

L  !  s-m  I  . 

k  =N 


L  =  K—N. 


PxM  = 


a  l  +  1 


a:  +  1 


x  >  0,  and 


W  ‘  vrrwf't',i  y  '  a 


p,(v)  =  j  *  Px(y  -  a)  Py(a)da 


- — - -  exp  -  f-f  V  ot1  1  exp  <  -  — — 1 —  [da.  (A7) 

(L  -  UKa;  +  1)  L  +  1  J  0  l  +  1  J 


Set  X  =  a/77.  Then 


Pji l)  = 


mat 


t]1  exp 


T3 -  M  X7  1  exp 

a:  +  1  0 


a;  +  1 


y<I\.  ( A8) 


